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Abstract 

It is proposed that superoscillations play an important role in the 
interferences which give rise to superluminal effects. To exemplify 
that, we consider a toy model which allows for a wave packet to travel, 
in zero time and negligible distortion a distance arbitrarily larger than 
the width of the wave packet. The peak is shown to result from a 
superoscillatory superposition at the tail. Similar reasoning applies to 
the dwell time. 
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g . 1 Introduction 

Superluminal effects have been predicted in conjunction with various quan- 
tum systems propagating in a forbidden zone. In these regimes, the proba- 
bility density current is imaginary, and various operational definitions of the 
velocity of a wave packet have been proposed - in many examples, giving 
differing values. In recent years, a number of experiments with superluminal 
photons have been performed, reviving interest in the problem, as well as 
controversy. The theoretical predictions have been verified - in fact, different 
tunneling times have been measured in accordance with the different op- 
erational definitions appropriate for the various experimental set-ups. (An 
excellent review on the experimental situation, as well as the theoretical 
background is afforded by 

Causality is not directly violated. The signal velocity (the velocity of 
propagation of an abrupt disturbance) is always subluminal. Other "veloc- 
ities" may well be superluminal, for example the group velocity of a wave 
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packet. In the latter case, the (local) peak of the packet appears at a point 
where constructive interference builds up, much earlier than the arrival of a 
freely propagating wave function. Thus the information stored in the tail may 
be travelling way ahead of the peak and can possibly be used to anticipate 
it. 

The purpose of this article is to investigate further the nature of the 
interferences that give rise to superluminality. It has been noted by Steinberg 
|fj, |3| that the superluminality phenomena is associated with post-selection: 
for instance, from a sample of particles that scatters of a barrier we examine 
only the subsample that tunnel through. When a pre- and post-selected 
system is subjected to a non-disturbing, "weak" measurement, the outcome 
of the measurement, known as the "weak value", can attain values that lie 
outside the spectrum of eigenvalues of the measured observable f§]. Weak 
values may hence be naturally related to the superluminal phenomena, as 
indeed, Steinberg has already argued that the dwell time is a weak value 
of a projector to the tunneling domain. The appearance of unusual weak 
values has been associated with a unique interference structure , for which 
Berry || coined the term " superoscillations" . As an instructive example of 
a superoscillatory function F(k) consider 



F(k; N, L) 



L/x \ jkx /N i ( 1 + L/x \ ikxo/N l h 



Here, N > 1 is an integer, and L and Xq being the super and reference shifts. 
For small k we expand exp(ikxo/N) and find: 



F(k;N,L) = e 



-ikL 



1+ ^ ^ +0(N->) 



e~ ikL . (2) 



Although this function is a superposition of waves e with \x\ < Xq, in the 



interval \k\ « / y L 2 — Xq = Ak, F{k) behaves nearly as a pure wave 
e %kL w j^h l arbitrary larger than Xq. In the regime \k\ < Ak the function 
oscillates rapidly. The number of these "superoscillations" is ~ \^N. This 
remarkable feature is derived at the expense of having the function grow 
exponentially in other regions. In the example above, for \k\ > Ak, we get 

In this article we will suggest that at least for certain cases, the construc- 
tive interference giving rise to superluminal effects, originates from a similar 
structure of superoscillations. To exemplify that, we consider a toy model, 
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(which extends on a previous proposal of Olkhovsky Recami and Salesi [0]), 
which allows for a wave packet to traverse, in a vanishing time and negligible 
distortion (the transmitted wave packet is the first derivative of the incom- 
ing packet), a distance arbitrarily larger than the original size of the wave 
packet. Hence the peak is here reconstructed from the exponentially small 
tail of the wave function. As far as we know, in the examples discussed to 
date, the superluminal shift of the wave packet is restricted. It is comparable 
to or much smaller than the initial wave packet size. We then show that the 
resulting null group-delay and dwell timesQ Finally, coming back to the role 
of post-selection we provide a rigorous proof for Steinberg's claim |fj, |3|] that 
the delay time is a weak value of a projector operator. 

The article proceeds as follows. In section 2 we calculate the dwell time 
and the group delay time for tunneling through n delta- function barriers. In 
the low energy limit, both turn out to be zero. We also derive the condition 
for the calculations to apply for a wave-packet. Using this condition, we see 
that for this system, the (negative) delay can be larger than the uncertainty 
associated with the length of the wave-packet. Section 3 deals with the rela- 
tion between superluminality and interference effects in the tail of the wave 
function, and superoscillations. The applicability to the example of section 
2, of the explanations of superluminality given in other cases, is discussed. 
Finally, in section 4 we elaborate on Steinberg's claim that the dwell time is 
a weak value. 

2 The group delay and dwell time for a par- 
ticle tunneling through an n-delta-function 
potential 

Olkhovsky Recami and Salesi [0] showed that a Schrodinger particle tunnel- 
ing through a double rectangular barrier traversed the distance between the 
bumps instantaneously in the limit that its kinetic energy was much smaller 

1 We follow the terminology of M- The group delay is defined g(i), where 

t is the transmission coefficient for tunneling through the region. Sometimes however, the 
same name is used for the difference between this value and the time it would take the 
same packet to traverse an equal distance in free propagation (i.e. for the additional delay 
introduced by the barrier). Other names for it are: phase- Wigner- and stationary phase- 
time. In the following we refer to the forward conditional dwell timc|^|, ||, as the dwell 
time. 
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than the height of the barrier. Unlike previous examples of superluminal tun- 
neling, the length of the region of superluminality consists of an arbitrarily 
long portion with zero potential, between the bumps. Replacing the rectan- 
gular barriers in the example discussed in by delta-function potentials, the 
calculations can be made somewhat simpler, and are easily generalized to n 
arbitrary delta bumps (still using the approximation of low kinetic energy). 

In this section we make a direct calculation of the transmission coefficient 
for the stationary scattering of a scalar particle obeying the Schrodinger equa- 
tion, off a multiple delta-function potential. The time independent equation 
is the same as for the scalar relativistic wave equation, and we focus on the 
Schrodinger equation merely for a simple concrete interpretation. 



2.1 Transmission through a multiple delta-function po- 
tential. 

Consider the Schrodinger equation with the following potential: 

V{x) = Eai5{x - Li); L = 0. (3) 
The energy eigenfunctions have the form (for x < and x > L n ): 



ip(x) 



Ae ikx + Be -ikx x < o 

Ce ikx + De~ ikx x > L n 



The coefficients satisfy: 




M=n 



i=l 



HkLi 



Pi 



1 e 
z — 1 



Pi 



mocj 
ik 



(4) 



(5) 



(6) 



In the limit of small kinetic energy (|/3j| ^> 1), we can drop the I matrices, 
as long as n < Pi. It is then straightforward to prove by induction on n that: 



M = ]JP t 



n? =2 (i - z t ) Yiuizi 1 - 1; 



, v -n?= 2 (i-^) -nu& -i; 

where z 1 — 1, Zi — ^KU-U-i) ^ _ i..n) 



0(1) 



(7) 
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As usual, we examine the case of "stationary scattering". To get the 
(amplitude) transmission coefficient for probability current flowing from the 
left,£, we put A = 1, B = r, C = t, D = in eq. (§: 

and we see that t = , so: 



A '"'fin o-l o-l Hi=2 Z ; 



U=2 



i=2 *i 



Tin i-\ ~ '\ "1 ^ n —I I 

n i=2 (i-^ n n =2(Zi 2 _ Zi2) 

J] 0-1 e -ifcZ« 

The stationary phase formula for the delay time, r g : 



(9) 



r ff = ^arg(i), (10) 

yields the value r g = — j£ = — for the delay, which cancels the time for 
a free particle, and we get an overall zero time for tunneling. Since this is 
true for all k, it should be true for an arbitrary wave packet, as long as the 
stationary phase approximation holds. The condition for that is derived in 
the next subsection. 



2.2 The condition for superluminal tunneling of a packet 

Restated for wave packets, our results so far can be summarized as: 

_ f / g(k) (e ikx + r(k)e- ikx ) e'^Hk x < ( . 

^' ' ~ { K-i^Ci^gi^e^-^e-^^dk x > L n ' ^ ' 

and 

TT /3 rl 

CAk\ = M2l 

1 ) -ik(-2i)"-i nr= 2 sin(A:(L, - L^)) 1 ) 

When Ak is sufficiently small, the diffusion can be ignored and C(k) 
can be considered constant (as will be shown shortly). Then we can again 
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separate out the time dependence of the wave function and the spatial part 
can be written: 



*V) = \ rJ {X) r ^ X< r (13) 
ry ' \ C(j)'(x - L n ) x > L n y J 

where <f>(x) in the two regions is related through analytic continuation. 
If 4>{x) = R(x)e lS ^ where R(x) is large and slowly varying in the region 
\x — x \ < Ax and S(x) goes through a few cycles there, then the time of 
arrival distribution of the transmitted packet will be approximately that of 
the incoming one, shifted by -F^r- Note also that this is also true for a mixed 
state which can be decomposed into various pure states with this property. 
Let us now find the explicit condition for C(k) to be approximately constant, 
for the case where Lj = ^L, aj = a, and as before, n < \(3\ = n f-. In this 
case, we have 

(-T 

— Ilk {2t sin kL/n) 
Using the fact that ^ = 1 + ^ + 0(x 4 ), we get: 

Thus, C(k) will be approximately constant if the spectrum of the wave packet 
is limited to k such that \k\ w In other words, A A; xs or Ax w 

^-j . This means that the length of the barrier can be arbitrarily longer than 
the "length" of the tunneling packet as usually defined (standard deviation 
of the x coordinate), the penalty paid being an exponential suppression of 
the amplitude. In passing notice the function F(k) in eqs. (1-2) displays a 
similar behavior. 



2.3 Calculation of the dwell time 

It is interesting to compare the "group delay" (which is zero in the low k 
limit) with the dwell time. For the sake of simplicity, we shall deal with the 
case n = 2. A direct calculation of the dwell time of the transmitted com- 
ponent can be made by calculating the transmission coefficient after adding 
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a potential which is constant over the region between the delta spikes, and 
vanishing outside it. We get: 

„— ikL 

t ~ /r 1 - , 

— 2% sin k'L ' 

where k! = ^2m(E — Vo)/h and Vo is the value of the potential between 

the deltas. Clearly, 9 = 0, and the (conditional) dwell time is zero as 
expected. 

3 Superluminality and its relation to inter- 
ference in the tail of the wavefunction 

The calculation of the transmission coefficient, t, can also be done in a way 
more suggestive of superoscillations. Let us explain this for the case of 2 
delta functions (the n = 2 case in ([|), "Fabry Perot interferometer"). 

Suppose a quasimonochromatic wave packet with wave number k arrives 
at the first delta spike. The transmitted component is the same as the origi- 
nal wave, except for an attenuation and phase which are independent of k. At 
the second delta spike, the wave splits into a (approximately unattenuated) 
reflected wave and a transmitted one which is apart from a fc-independent 
multiplicative constant the same as the impinging wave. The reflected com- 
ponent is again reflected at the first delta, and arrives at the second delta 
with an additional phase of 2kL, but with approximately the same ampli- 
tude as the original transmitted wave. In a like manner, one gets additional 
transmitted waves with additional phases of 2nkL,n = 2,3,..., and ampli- 
tudes which decrease very slowly. Thus we get the following formal sum for 
the resulting amplitude of the wave (up to a multiplicative constant): 

E ikxi ik1L\n _ _ ikx 

a 1 ' -2z sin AlL 1 ' 

which is in agreement with our previous calculation. This is an example of 
superoscillations since a sum of positive wave vectors rezults in a negative one 
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(or, equivalently, a sum of positive shifts which results in a negative shift). [] 
This is true in the following sense: for \k\ << 1/L the denominator of the 
right hand side can be considered constant. However, in such a small interval 
the function does not really oscillate, so it really doesn't have a well defined 
frequency. To really speak about superoscillations we need to have a large 
number of delta-functions. The sum for the case n > 2 factors into n— 1 sums 
of the above form, in the low kinetic energy limit, since to leading order in 
/3 , the only contributions are from waves which tunnel through each delta 
only once, but may be reflected any number of times between consecutive 
deltas. We then reproduce the results of subsection 2.2, where we had the 
weaker condition k < \/n — 1/L, allowing the function to complete many 
oscillations in the region. Note the similarity to the situation described by 
eqs. (1,2). 

A different calculation of the (conditional) dwell time than that of section 
2.3 assumes a "clock" which is activated by the presence of the particle in 
the region [0, L], i.e. a degree of freedom r with an interaction Hamiltonian 
-f^int = @{o,l](x)Pt, where p T is the canonical momentum conjugate to r. If the 
initial state of the clock has small enough Ap T , one gets for the final phase of 



r the same expression as eq. I6|, with x, k replaced by the clock's coordinates. 
This can be interpreted as the sum appropriate for a weak measurement of 
@[o,l]{%), as shown in the next section . Note that in this case, not only do 
the dwell and phase times coincide, but they are also described by the same 
mechanism. 

The group delay in tunneling through a thick barrier follows from the 
fact that under the barrier, no phase accumulates, and the entire phase shift 



2 The sum in eq.(16) actually diverges, the physical reason being that we have neglected 
the attenuation of the amplitude, in order to maintain consistency with the low kinetic 
energy approximation we have used so far. For the case n = 2 it is easy to evaluate eq.(^|) 
without resort to that approximation, and the resulting transmission amplitude is: 

0- 2 

m = 7 7^77 ( 17 ) 







_ e 2ikL 









Similarly, the sum on the left of eq. dl6|) should be replaced by the exact one: 



(i+/r 2 £ 

3=0 L 







1 + /3 



2 V 

2ikL 



0~ 2 

(18) 





- 1 - 




p2ikL 
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comes from the boundaries and is practically independent of the thickness. 



For cases where interference with a delayed wave, a few authors [12, 13, 14, 15 



have suggested a different explanation. In Chiao and Steinberg's words |L|: 
"If destructive interference is set up between part of the wave travelling 
unimpeded and part which has suffered a delay At due to multiple reflections, 
one has ^ out (t) = * in (t) At) « (1 - £)% n (t) + £Atd% n (t)/dt « 

(1 — £)^i n (t + £At/(l — £)), which is already a linear extrapolation into the 
future. In cases where the dispersion is sufficiently flat, as in a bandgap 
medium, the extrapolation is in fact surprisingly better than this first-order 
approximation. As was suggested by Steinberg |b|] and recently discussed 



more rigourously by Lee and Lee [16] and Lee [17], this implies that even a 



simple Fabry-Perot interferometer exhibits superluminality when excited off 
resonance" [presumably, £ << At ^,^ )■ We would like to explain this "better 
than first-order" approximation. Let us instead look at the momentum wave 
function. A spatial shift corresponds to a linear shift in this function. A 
positive spatial delay would correspond to a linear shift steeper than one, and 
the converse for a negative delay. In the Taylor expansion of the transmission 
coefficient for the momentum wave function, the zero term is insignificant, 
the second corresponds, as just explained, to the spatial shift, and the higher 
give the distortion. When many large and evenly distributed shifts waves 
interfere, their sum is for a wide range of momenta, zero, and in particular 
momentum independent. In other words, the momentum wave function is 
flat for a wide band of frequencies. This corresponds to a much better than 
first order approximation of the spatial wave function, as can be seem in the 
special case of the system of section 3 of this paper. 



4 The dwell time as a weak value 

We would like to calculate the expectation value of the time measured by 
a "clock" consisting of an auxiliary system which interacts weakly with our 
particle as long as it stays in a given region. Furthermore, we would like to 
restrict the calculation only to the subensemble of particles which ultimately 
end up on the right of the barrier. The simplest interaction is perhaps the 
one defined by the Hamiltonian: 

H int = P r X (0,L) 
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Where r is the degree of freedom of the clock and P T is its conjugate mo- 
mentum, and X( 0j l)(x) = j q ^otherwise^ ' * S ^ e e ^ ec ^ ve f° rm 5 
for example of the potential , seen by a particle in an S z eigenstate, in the 
Stern-Gerlach experiment (r being the z coordinate, and (0, L) the region 
of the magnetic field). Assuming the clock to have at large negative times 
the expectation value 0, a perturbation calculation shows that following the 
interaction with the particle and the subsequent post-selection of the particle 
state to be |/), the expectation value of r at large positive times is given by 
the formula: 

, N I-oo dt Io ' dxW f (x,t)Vi(x,t) . 
E(r,t^oo|i,/) = C(M . (Xi0) » j(Xi0) (19) 

Steinberg^, |3| has arrived at this formula under similar assumptions by 
a somewhat different line of reasoning. He introduced the term conditional 
(quantum) probability for the probability distribution of a system following 
post-selection, and we use the notation of conditional expectation in the 
formula above, in the same spirit. As noted by Steinberg, the last equation 
is a special case of a weak value. 

This formula is valid when r and p T do not appear in additional terms 
in the full Hamiltonian, but the generalization is straightforward. To prove 
the formula, let us work in the interaction picture. Denote the initial state 
of r by \(p T (t)) (and the corresponding wave function by <p(r,t)), and the 
initial (preselected) and final (post-selected) states of the tunneling particle 
by \i(t)), \ f(t)) (^i(x,t), ^f(x,t)), respectively. Using the interaction picture 
and expanding to first order in 

\<j>r,t)\i(t)) = Te-*f-°° Hl{t ' )dt '\<i) T} t^ -oo}\t(t -> -oo)) ~ 



i 



'-a j. 



f H^dA |0 T ,*-> -oo)|i(t -> -oo)) (20) 



where 

^mt 



3 We can satisfy the condition of small perturbation by choosing an initial clock wave 
function concentrated about P T = with small enough uncertainty. 
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After making the post-selection of state |/) for the particle, the clock is left in 
the state given by the above expression, multiplied on the left by ^~^~~a ■ 
The expression we get after writing out the explicit form of Hi nt is: 

<f>(r,t -y +00) ~ e~ l{x (°^ )wPT/h (j)(T,t -> -00) 

= <f>(r + (X {0tL) ) w ,t ^ -oo) (21) 

where 

{X{0 ' L))W ~ C« K x,0)^(x,0) (22) 

(the integral in the denominator is evaluated at t = merely for convenience 
- it is of course, time invariant). The wave functions in the integral in 
the numerator can be taken in the Schrodinger representation. And the 
expectation for r at large positive times is shifted (from its initial value of 
0) by this value, as claimed. In contrast to the familiar eigenvalue spectrum 
of a physical operator, its weak values can take any complex values.^ 



The research was supported in part by grant 62/01-1 of the Israel Science 
Foundation, established by the Israel Academy of Sciences and Humanities, 
NSF grant PHY-9971005, and ONR grant N00014-00-0383. 

4 To see this, let us develop the initial and final states of the particle in terms of the 
eigenfunctions of the operator to be measured, A: 

\i >=^2a k \a k >,\f >=^2(3 k \a k >, (A\a k ) = a k \a k )) (23) 

k k 

Then, Aw — ^fpr = ^" fcafc If A is nontrivial, it has more than one eigenvalue. 
Assume k — 1,2 correspond to two of these, and take Pi = P2 — ^ Then the two 

equations: Aw — Ql °^°^ a2 = z, and |ai| 2 + |a2| 2 = 1, are 3 real equations in 4 unknowns. 
They can be solved for any value of z, as can be verified easily. 
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